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Summary. Globular clusters are often assumed to be good tracers of major star 
formation episodes in their host galaxies. While observations over the past 2 decades 
have confirmed the presence of young objects with globular cluster-like properties in 
many galaxies, it is still not well understood exactly how the formation efficiency 
of bound star clusters relative to field stars and their mass spectrum depend on 
external factors. The cluster initial mass function typically appears to be consistent 
with a power-law with a slope ~ —2, but most attempts to constrain any upper 
limit on the CIMF have been limited by size-of-sample effects. However, evidence 
is starting to accumulate for possible truncation of the cluster mass function. It is 
tentatively suggested that the upper mass limit may currently be at ~ 10^ M© in 
the Milky Way disk, while there are indications that it is ~ 5 x 10"" Mq in M51 
and about 2 x 10^ Mq in the Antennae. Some extreme starbursts (e.g. Arp 220, 
NGC 7252) are (or were) able to form clusters as massive as ~ 10'' M0. The overall 
formation efficiency of star clusters (relative to field stars) in the Galactic bulge 
may not have been much different from that in the disk today, but was probably 
significantly higher for metal-poor GCs in halos. 

1 YMCs — Guides to Young Stellar Populations? 

Three or four decades ago, globular cluster (GC) formation was thought by 
many to be a phenomenon occurring only in the early Universe (e.g. |2U|). In 
the meantime, young, compact star clusters with masses in the range 10''- 
10^ M0 have been found in many different galaxies (e.g. |18[ I21L 1^ ). and 
there is a growing concensus that these "Young Massive Clusters" (YMCs) 
may well be young analogues of the old GCs associated almost universally 
with the spheroidal components of galaxies. By implication, GCs have become 
potentially interesting not just as fossil left-overs from the early Universe, but 
more generally as test particles for studies of extragalactic stellar populations. 

In keeping with this spirit, contemporary observing proposals or papers on 
GCs often include an introductory remark along the lines of: "GCs are thought 
to be good tracers of the major star forming episodes in their host galaxies ". 
YMCs/GCs may indeed trace star formation quite generally, but it is also 
clear that some caution must be exercised. For example, the number of GCs 
per field star varies from galaxy to galaxy (the classical "specific frequency 
problem"), as well as between stellar populations within galaxies ^lE]. This 
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must be due to differences in the formation efficiency of GCs relative to field 
stars, in the GC survival rates, or (more likely) some combination of the two. 
Of the roughly 150 catalogued GCs in our own Galaxy about 2/3 are 
associated chemically, spatially and kinematically with the halo, while this 
is true for only ^ 1% of the stars. Conversely, some 90% of the stellar mass 
resides in the thin disk of our Galaxy, while no GCs are currently known to 
be associated with this component. Even though some quite massive clusters 
might be located in remote parts of the Galactic disk |3 EH , a simple scaling 
by stellar mass of the GCs in the bulge or halo would predict hundreds or 
thousands of GCs in the disk, seemingly at odds with the observations. 

All this is of course well known, and is the reason why GCs are often 
assumed to trace only "major star forming episodes". The problem then re- 
mains to define when a star forming episode qualifies as "major". Perhaps 
GCs mainly trace spheroid formation ^ , but some stars currently residing in 
spheroids may originally have formed in disks. As an example, it is illustra- 
tive to consider the outcome of merging two Milky Way galaxies: The merger 
product would contain about 200 metal-poor and 100 metal-rich pre-existing 
GCs from the progenitor galaxies, assuming no GCs are destroyed. The cur- 
rent gas mass in the Milky Way is about 0.5 x 10^° Mq [Hj, which is about half 
the mass of the bulge. Assuming that the merger would form GCs with the 
same efficiency as the bulge, using all the available gas, about 50 new metal- 
rich GCs would form (and survive). The resulting GC population would then 
consist of three major sub-components: metal-poor, old GCs originating in 
the progenitor galaxy halos, moderately metal-rich GCs from the pre-existing 
bulges, and very metal-rich GCs formed in the merger. The estimate of the 
number of new GCs is obviously very crude. However, the main point here 
is that while the majority of the stars in the resulting spheroid (about 90%) 
would have formed in the disks already before the merger took place, their 
formation history might not be reflected in the GC system. 

Mergers at higher z were likely more gas-rich, and the discrepancy between 
the star- and GC formation histories may have been less extreme than in the 
example outlined above. Nevertheless, since GCs play such an important role 
in attempts to constrain the star formation histories of early- type galaxies, 
there is a clear need to also quantify the limitations better. 

2 What can we learn from studies of YMCs? 

In order to understand the differences between properties of star cluster pop- 
ulations in different galaxies better, it is useful to divide the problem into 
three sub-problems which can, to a large extent, be addressed separately: 1) 
Understanding the cluster (initial) mass function (CIMF) - Is it universal, or 
do some parameters (e.g. the slope or upper mass limit) vary as a function of 
external parameters (star formation rate, gas pressure/density)? 2) The clus- 
ter formation efficiency relative to "field" stars: Again, how does this depend 



Young Massive Clusters 3 



on external factors? 3) Disruption effects: these are driven both by internal 
mechanisms (two-body relaxation, binaries, black holes) and external factors 
(shocks, interactions with giant molecular clouds, tidal fields). However, it is 
a hard problem and progress has been slow in making reliable, quantitative 
predictions for the time scales involved, at least until very recently. 

In the following I will concentrate mainly on the first item in the list, 
with only a few remarks about formation efficiencies. Disruption effects are 
covered elsewhere in this volume (e.g. Baumgardt, De Marchi, Vesperini). 

2.1 The Cluster Initial Mass Function 

The number of galaxies with reliable constraints on the CIMF remains small. 
Probably the best-known example is the Antennae, where the CIMF appears 
well approximated by a power-law dN/dM oc with a « — 2 over the 
range 10'' < M/Mq < 10^ ^7\. Similar mass functions have been found in 
M51 0, NGC 3310 and NGC 6745 13, the Milky Way 9 and in the LMC 
|15|. although not all studies cover the same mass range. It seems reasonable 
to conclude that star clusters typically form with a mass spectrum that can be 
well approximated by a uniform power-law over some mass range. It should be 
mentioned that some dwarf galaxies contain a few clusters which are much 
brighter than one would expect from the total number of star clusters in 
those galaxies |2j. Here, however, I mainly focus on the opposite problem, 
i.e. whether there might be a truncation of the CIMF at some upper mass 
■^trunc, thus inhibiting efficient formation of clusters above a certain mass 
limit (Mtrunc) under certain circumstances (e.g. in the Milky Way disk today). 

2.2 Some considerations on the Milky Way 

Starting again with the Milky Way, it is interesting to consider the conse- 
quences of postulating that young clusters are drawn purely at random from 
a power-law distribution with a = —2. The current star formation rate in 
bound star clusters in the solar neighbourhood is estimated to be around 
5.2 X 10~^° M0 yr~^ pc~^ ^B]- Assuming for simplicity a constant cluster 
formation rate in the Galactic disk over the past 10 Gyrs and that the cluster 
formation rate in the Solar neighbourhood is representative for the disk as 
a whole, this corresponds to ~ 10^ M0 formed in bound clusters within the 
Solar circle (8.5 kpc). This is most likely a conservative estimate and the 
true number may well be significantly higher. Sampling these clusters from 
a power-law CIMF for 10^ < M/Mq < lO"^, one predicts a total of close 
to 9 X 10^ clusters formed over the lifetime of the disk, of which 9000, 900 
and 80 have M > IO^Mq, M > IO^Mq and M > IO^Mq. These numbers 
do not depend strongly on the adopted upper and lower integration limits, 
although they do depend on the CIMF slope - a steeper slope implies a more 
bottom-heavy CIMF, with fewer high-mass clusters. 
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There are two consistency checks worth making: first, the steUar mass of 
the Milky Way bulge is about 10% of that of the disk. If clusters formed 
with the same efficiency in the bulge as they do in the disk now, one might 
expect about 90 clusters with M > IQ^Mq in the bulge, of which 8 have 
M > IO^Mq. The actual observed numbers of GCs in the bulge are smaller 
by a factor of 4, which may be partly due to disruption effects. However, there 
is no indication that clusters formed with a higher efficiency in the bulge, and 
the number of GCs observed in the bulge may be roughly consistent with a 
formation efficiency (relative to field stars) similar to that observed today in 
the Galactic disk. The numbers for the halo are more difficult to reconcile 
with this picture, since the halo has about twice as many GCs as the bulge 
but an order of magnitude fewer stars. This suggests a higher formation 
efficiency of metal-poor halo GCs, consistent with observations of early-type 
galaxies |14[ The distinction between halo (metal-poor) GCs and all 
other star clusters may be more fundamental than the one between old GCs 
in general and present-day star formation. Note that these arguments are 
slightly different from those of McLaughlin ^HI who argued for a universal 
cluster formation efficiency relative to the total available gas mass. 

Second, we can compare with the number of massive star clusters actu- 
ally observed in the disk. By construction, the current formation rate agrees 
with the number of low-mass clusters {M < 10^ M©) observed locally, but 
it is interesting to see what happens when extrapolating to higher masses. 
Recently, at least two young clusters with masses in the range lO^'-lO^ M© 
and ages < 10'' years have been identified |Z|EI|. If they are taken as rep- 
resentative of the formation rate of such objects within a distance of 5 kpc, 
this corresponds to 2 kpc~^ Gyr~^, or 4500 such clusters formed within the 
Solar circle over 10 Gyrs. Again, this is roughly consistent with the order-of- 
magnitude estimates above (in fact, 8000 clusters with 10"* < M/Mq < 10^ 
are predicted), suggesting that clusters like Wester lund 1 with masses up to 
about 10^ Mq occur naturally (albeit rarely) as part of the normal hierarchy 
of star formation in the Milky Way disk today. 

A typical disruption time for a cluster with mass M in the solar neigh- 
bourhood is tdis = 1.3 Gyrs (M/10*Mo)° '^2 ^ Assuming that this scaling 
remains valid for M > 10'* M©, a cluster with an initial mass of 10^ M© 
is expected to disrupt in about 5 Gyrs while a 10^ M© cluster has a life- 
time well in excess of a Hubble time. The disk should then contain about 
500 clusters with masses greater than 10^ M© and still virtually all of the 
80 clusters with M > 10^ M© formed over its lifetime. Of these objects, 7 
should be found within a distance of 1 kpc. This estimate makes the crude as- 
sumption that clusters disrupt instantaneously, while in practice mass is lost 
at a nearly constant rate over the lifetime of the cluster. Although current 
catalogs of Milky Way open clusters are highly incomplete beyond 1 kpc, it 
seems unlikely that a large population of clusters with masses in the range 
10^-10^ M© could have been missed. It appears plausible that the CIMF 
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in the Milky Way is truncated somewhere in the vicinity of 10^ Mq, or at 
least becomes much steeper than a power-law with a slope of —2. However, 
it would be highly desirable to quantify better the completeness of current 
cluster surveys in the disk out to large distances (several kpc). 

2.3 Constraints on the CIMF in other galaxies 

Studies of the CIMF in external galaxies are complicated by the rapid change 
in mass-to-light ratio that characterizes simple stellar populations. Observed 
luminosities cannot be converted to masses unless the age of each individual 
cluster is known with some accuracy, which generally requires U-band imag- 
ing. Several studies have shown that the luminosity function of young star 
clusters generally appears to be sampled all the way up to its statistical upper 
limit 1^ 121 El ■ Even the scatter around the predicted relation is consistent 
with random sampling |17) . No direct evidence for truncation of the LF has 
been found so far, suggesting that most large galaxies are physically able to 
form star clusters with masses at least up to about 10^ Mq. 
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Fig. 1. Mass and My distribution for a simulated cluster population with 3 x 
10^ clusters and 6<log(age/yr)<9. Hashed and outlined histograms are for mass 
distributions truncated at 10^ Mq and 10* Mq. 

In general, there is no straight-forward way to infer the mass function of a 
cluster sample from the observed luminosity function 27 . Only in the special 
case where the MF is a simple, untruncated power-law, or the age distribution 
is a delta function, will the LF and MF have the same shape. This point is 
illustrated in Fig. ^ which compares the mass- and luminosity functions 
for simulated cluster samples with power-law mass functions truncated at 
10^ Mq and 10^ Mq. The clusters were assigned random ages uniformly 
distributed between 10^ years and 10^ years and masses were converted to 
My magnitudes using SSP models 0. Mass loss and cluster disruption were 
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ignored. In the left-hand panel, the MF truncation occurs as a simple cut- 
off, while in the LF (right) the mass cut results in a steepening of the slope 
at the bright end rather than any distinct cut-off in luminosity. For a real 
cluster population, the maximum mass may be inferred from the "bend" that 
occurs at My ~ —8 in Fig. assuming that the MF is a power-law up to 
-Mtrunc Such a bcud has been observed in the Antennae and M51, where 
it may be explained by a MF truncated near 2 x 10^ Mq and ^ 5 x 10^ Mq 
123 El- However, in still more active galaxies such as Arp 220, NGC 7252 
and NGC 1316, there are clusters with masses as high as ~ 10'' Mq [T1 I26|. 

Let us now consider the behaviour of the following quantities as a function 
of the total number of clusters (N) in a population with a truncation at an 
upper mass limit Mtiunc: 1) the maximum mass Mmax occurring in the popu- 
lation, 2) the mass of the brightest cluster Afbrightcst, and 3) the magnitude of 
the brightest cluster, j\/^rightcst ^ ^ clusters are sampled at random from a 
power-law dN/dM oc with lower mass limit Mmin, then statistically the 
most massive cluster will have Mmax = N M^ain jMHH- From Monte-Carlo 
simulations of various cluster formation histories, Weidner et al. found 
that when clusters are sampled at random from a power-law MF, the most 
massive cluster is also the brightest in about 95% of the cases. If the mass 
function has a real physical upper limit, this is not necessarily the case. 

The left panel in Fig. |21 shows the results of Monte-Carlo simulations for 
-^max and Afbrightcst in cluster samples with various Mtrunc limits. Clusters 
were drawn at random from the same population used in Fig. ^ truncated 
at Mtrunc values between 10'* Mq and 10^ Mq. The median values of Mj^ax 
and Afbrightcst in 1000 experiments are shown with solid and dashed lines 
as a function of the number of clusters sampled. The relation predicted by 
pure sampling statistics is indicated as a dotted line. For small TV and large 
A/trunc, the truncation of the MF is not "felt" and the Afmax and Afbrightcst 
vs. iV relations approach the curve for random sampling, i.e. the Weidner et 
al. result is reproduced. Conversely, for small Aftrunc and large N, statistical 
effects become unimportant and Afbrightcst ~ Af,„ax ^ Aftrunc- In between 
these extremes there is a regime where the mass function is likely to be sam- 
pled up to higher masses than the luminosity function, or Afmax > Afbrightcst- 
This turns out to be the situation encountered in many real galaxies. 

For both Afmax and Afbrightcst, Fig- 121 shows a fairly rapid transition 
between the regimes where size-of-sample effects and truncation dominate. 
Thus, it might seem that the observed strong relation between total number 
of clusters in galaxies and the luminosity of the brightest cluster |21 EI 
is incompatible with truncation of the CIMF playing any important role. 
However, as shown in the right-hand panel of Fig. |21 j^g^g ^ gteep 

dependency on N over a much wider range in N than Afmax and Af brigiitcst - 
This is because the mean age of the brightest cluster shifts towards younger 
ages for higher N, as it becomes increasingly likely to encounter a cluster 
with mass near Aftrunc in the brief phase where the M/L ratio is very low. 
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Fig. 2. Left: Mass of tlie most massive (solid lines) and most luminous (dashed 
lines) cluster as a function of total number of clusters in simulated cluster sam- 
ples with various fixed upper mass limits. The dotted line is the maximum cluster 
mass expected from random sampling from an untruncated power-law. Right: The 
luminosity of the brightest cluster as a function of and Mtrunc- 



From the preceding discussion it should be clear that the LF can be domi- 
nated by sampling effects, even if the MF does have a physical truncation. For 
example, for Aftiunc = 10^ Mq, the Mmax curve starts to flatten at ~ 10^, 
implying that the MF would be sampled up to its physical upper limit al- 
ready in a relatively cluster-poor galaxy. The ]\,fy"^^^°^^ curve, on the other 
hand, would continue to rise beyond N = 10^, corresponding roughly to a 
Milky Way-sized galaxy (cf. Section . 

The key to detecting a physical upper limit of the CIMF will be to cover 
a dynamic range large enough to study the LF in detail and detect signatures 
such as the "bend" in Fig. ^ Ultimately, however, inferences about the MF 
based on the LF remain indirect and dependent on assumptions about the 
shape of the MF, star formation history etc., and ideally it would be desirable 
to have direct information about the mass functions in several more galaxies. 
Several such studies are now underway, and may provide important insight 
into the MF in the not too distant future. 



3 Concluding remarks 

While research in extragalactic star clusters remains a very active field, we 
are still facing a number of important questions. The ubiquity of YMCs in 
external galaxies is now well established, but the apparent absence of a large 
population of clusters with masses in the range 10^-10^ Mq in the Galactic 
disk remains somewhat of a puzzle. This may suggest an upper limit to the 
CIMF near 10^ Mq in the Milky Way, as noted already by van den Bergh 
& Lafontaine [^SIj although studies of disk clusters are still hampered by 
our own location close to the Galactic plane. However, the recent realization 



8 S. S. Larsen 



that the Milky Way is stiU forming clusters with masses near 10^ Mq hints 
that the CIMF in the Milky Way may not be very different from that in the 
LMC, the main difference being that the census of YMCs is more complete 
in the LMC. There is some evidence for truncation at a somewhat higher 
mass in M51 5 X 10^ Mq) and the Antennae (~ 2 x 10^ Mq), while Arp 
220, NGC 1316 and NGC 7252 host clusters as massive as - lO'^ Mg. These 
galaxies also define a sequence of increasing star formation rate, suggesting 
that galaxies with higher SFRs are physically able to form more massive 
clusters. This may provide a hint as to why GC formation was common at 
high z, when SFRs and gas densities were generally higher. 

References 

1. N. Bastian, R.P. Saglia, P. Goudfrooij, et al.: A&A, 448, 881 (2006) 

2. O.H. Billet, D.A. Hunter & B.G. Elmegreen: AJ, 123, 1454 (2002) 

3. A. Bik, H.J.G.L.M. Lamers, N. Bastian et al.: A&A, 397, 473 (2003) 

4. J. P. Brodie, J. Strader: Ann. Rev. Astron. Astrophys., in press (2006) 

5. G.A. Bruzual, S. Chariot: MNRAS, 344, 1000 (2003) 

6. B.W. Carroll, D.A. Ostlie: An Introduction to Modern Astrophysics, Addison- 
Wesley Publishing Company (1996) 

7. J.S. Clark, I. Negueruela, P.A. Crowther, S.P. Goodwin: A&A, 434, 949 (2005) 

8. R. de Grijs, P. Anders, N. Bastian, et al.: MNRAS, 343, 1285 (2003) 

9. B.G. Elmegreen & Yu.N. Efremov: ApJ, 480, 235 (1997) 

10. D.F. Figer, J. MacKenty, M. Robberto, et al.: ApJ, in press (,astro-ph/0602146j 
(2006) 

11. J.C. Forte, F. Faifer, D. Geisler: MNRAS, 357, 56 (2005) 

12. M. Gieles, S.S. Larsen, N. Bastian, & I.T. Stein: A&A, 450, 129 (2006) 

13. W.E. Harris: AJ, 112, 1487 (1996) 

14. W.E. Harris, G.L.H. Harris: AJ, 123, 3108 (2002) 

15. D.A. Hunter, Elmegreen, B.G., Dupuy, T.J., & Mortonson, M.: AJ, 126, 1836 
(2003) 

16. H.J.G.L.M. Lamers, M. Gieles, N. Bastian, et al.: A&A, 441, 117 (2005) 

17. S.S. Larsen: AJ, 124, 1393 (2002) 

18. S.S. Larsen: in: Planets to Cosmology: Essential Science in Hubble's Final 
Years, ed. M. Livio, STScI (2005) 

19. D.E. McLaughlin: AJ, 117, 2398 (1999) 

20. P.J.E. Peebles & R.H. Dicke: ApJ, 154, 891 (1968) 

21. F. Schweizer: In: "Extragalactic Star Clusters", ASP, eds. D. Geisler, E.K. 
Grebel, D. Minniti (2002) 

22. C. Weidner, P. Kroupa, S.S. Larsen: MNRAS, 350, 1503 (2004) 

23. B.C. Whitmore, Q. Zhang, C. Leitherer, et al.: AJ, 118, 1551 (1999) 

24. B.C. Whitmore: In: "A Decade of HST Science", eds. M. Livio, K. Noll, M. 
Stiavelli, Cambridge: Cambridge University Press, p. 153 (2003) 

25. S. van den Bergh & A. Lafontaine: AJ, 89, 1822 (1984) 

26. CD. Wilson, W.E. Harris, R. Longden, N.Z. ScoviUe: ApJ, 641, 763 (2006) 

27. Q. Zhang & S.M. Fall: ApJ, 527, L81 (1999) 



